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Abstract 

We suggest a method of summing the perturbation theory (PT) asymptotic 
series related to infrared (IR) renormalons in QCD by using special renormalization 
schemes in which the running coupling can be integrated over the region of small 
momenta. For the method to work, one should consider higher order PT corrections 
to the standard bubble-chain diagrams. High-order corrections allow one to choose 
a scheme in which the coupling evolution can smoothly be extrapolated to small 
momenta. In such schemes the sum of an (extended) IR-renormalon asymptotic 
series is defined as an integral of the running coupling over the IR region. We 
give explicit examples of renormalization schemes in QCD which can be used for 
summing IR-renormalon asymptotic series according to our definition. 



High-order behavior of perturbation theory series in quantum field theory determines 
the analytic structure of amplitudes as functions of the coupling constant. This analytic 
structure encodes important information on the ultimate precision obtainable in the finite- 
order PT analysis for physical observables. A general approach to analyzing high-order 
PT expansions is based on the saddle-point approximation for the path integral evalua- 
tion Jl], H, ||. The leading asymptotic behavior of the PT coefficients has been investigated 
in some models of quantum field theory while there are no exact results in QCD. The 
high-order behavior can also be studied by selecting large contributions of some particular 
subclass of diagrams 0, 0. In QED a set of diagrams generated by inserting an infinite 
number of one-loop fermionic corrections into the photon propagator ( "bubble-chain" di- 
agrams) is usually considered as a relevant one in this context. In other gauge models 
(for instance, in QCD) the selection of relevant diagrams is more involved due to the ne- 
cessity to retain gauge invariance. The divergent series |^ emerging from such a type of 
approximation for the high-order behavior - "bubble chains" or renormalons - have inten- 
sively been studied 0, ||, . The ultraviolet (UV) renormalons refer to asymptotic series 
emerging from the integration over large energies while the infrared (IR) renormalons cor- 
respond to integration over small energies. Depending on the model, renormalons lead to 
alternating (Borel summable) asymptotic series or non-alternating (Borel non-summable) 
asymptotic series. It is generally assumed that the Borel summation procedure gives the 
proper value for the asymptotic series which means that the theory can cope with the 
renormalon singularities that give alternating asymptotic series. The appearance of Borel 
non-summable series is usually considered as a signal that a theory (or an approximation 
in which it is studied) may be inconsistent. In QCD the IR renormalons give an example 
of Borel non-summable series if the leading-order approximation for the QCD /5-function 
with a positive coefficient Pq (not too large number of fiavors) is used. 

The recent interest in IR renormalons in QCD has been connected with nonpertur- 
bative (nonPT) power corrections This interest is caused by a limiting accuracy of 
finite-order PT expressions that cannot presently compete with a high precision of experi- 
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mental data (e.g. ||TT1); the accuracy of theoretical formulae has to be improved. A regular 
way of going beyond a finite-order PT approximation consists in the inclusion of nonPT 
terms as power corrections to the ordinary PT expressions. For the two-point correlators 
of interpolating hadronic currents the nonPT terms can be represented phenomenologi- 
cally as vacuum expectation values of some local operators |12|, |l^ within the operator 
product expansion (OPE) which is a general method of quantum field theory [jl4|, |l^, |16 



Such an approach has been successful in investigating hadron properties in QCD |]T] 



However, there is no well-established technique to calculate nonPT contributions in cases 
when OPE is not applicable. It has been suggested to consider the IR renormalons in 
QCD (or, generally, Borel non-summable asymptotic series of perturbation theory) as a 
source for power corrections in such cases. 

In QED, the bubble-chain diagrams form a gauge invariant class and give dominant 



contributions in the framework of 1/Np expansion [|l8l. In QCD the relevant diagrams are 
the gluon bubbles which do not form a gauge invariant set. Therefore, to find renormalon 
contributions in QCD one applies a method of "naive nonabelianization" to the same set 
of diagrams as in QED fl^, This is equivalent to the substitution of the leading or- 
der running coupling constant into integrals over the phase space of the original diagram. 
Thus, the renormalon approximations are attractive for discussing power corrections from 
a practical point of view because the relevant diagrams are closely connected with the 
running of the coupling constant within the renormalization group analysis. The renor- 
malon contributions can be obtained by either a direct summation of the selected set of 
diagrams or performing the integration of the running coupling constant over the phase 
space of the leading order diagram which serves as a source for a bubble-chain generation. 
These ways of obtaining the renormalon contributions (the direct summation of a given 
subclass of diagrams vs the integration of the running coupling constant) are similar to 
the methods by which the ultraviolet asymptotics of Green's functions in quantum elec- 
trodynamics was originally found. First the direct summation of the leading ultraviolet 
asymptotics in all orders was performed and later the same results were obtained by 
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computing the one-loop /3-function and solving the renormalization group (RG) equation 



for the invariant charge p2| , p3| . Note that to use the direct summation in high orders of 
PT is technically difficult while to solve the RG equations with a high-order /5-function 
is feasible. This observation is also valid for renormalons: to analyze corrections to the 
leading order renormalon expressions within the direct summation approach is technically 
difficult while the integration of the running coupling constant can practically be done at 
any order of the approximation for the evolution if the corresponding /3-function is known. 
The technique for calculation of RG functions is well developed p4| and the high-order 
PT expressions for /3-functions are available for a number of models. 

In the asymptotically free QCD, the running coupling is regular at large momenta 
and integrals over large momenta are well defined which solves the problem of summing 
the UV-renormalon contributions. In the direct summation technique these contributions 
correspond to alternating asymptotic series which are Borel summable and can be obtained 
by a formal term-by-term integration of the leading order PT expansion for the running 
coupling. For determining the IR-renormalon contributions in QCD one should integrate 
the running coupling over small energies where it is singular in the leading order of PT 
(Landau pole). Therefore, the integral (and the contribution of IR renormalons) is not 
defined in the leading order of the coupling evolution. By expanding the running coupling 
in PT series under the integration sign and term-by-term integration of the resulting series 
one arrives at a divergent series which is not Borel summable. This corresponds to the 
direct summation method. Thus, the problem of summing the contributions stemming 
from the infrared region can be reduced (at least at the formal computational level) to 
defining an integral of the running coupling over the small momenta. The interpretation 
of the integral is necessary if there is a singularity in the IR region as in the leading order 
of PT in QCD where there is a Landau pole. A particular definition of the evolution in 
the IR region provides a recipe for resummation of the Borel non-summable asymptotic 
series emerging from the integration over the small momenta. Because the evolution in 
the IR region is not perturbative the extrapolation of the running to small momenta is not 
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unique and implies some nonPT assumptions. For instance, the Landau pole appears for 
the evolution generated by the leading order PT /3-function. This singularity is unstable 
against inclusion of higher order PT corrections which are large in the IR region because it 
is a domain of strong coupling. Higher order corrections can strongly change the behavior 
of the running coupling at small energies. For instance, the pole can disappear (infrared 
fixed point) or the singularity can become an integrable one and the integration will lead 
to a physically acceptable result. In such cases the integration is possible that defines 
the way of resummation of the asymptotic series generated by the IR renormalons which 
means that high-order PT corrections stabilize the theory. An acceptance of a particular 
scenario for the IR evolution, however, cannot be justified in PT framework only and 
means that some nonPT assumptions have been made. 

In the present paper we define the renormalon contributions using the method of the 
integration of the running coupling. The relevant representation can be written in a 
compact form and generalized to higher orders of PT. In the leading order it is equivalent 
to the "naive nonabelianization" trick. The IR contribution in QCD is not defined in 
the leading order because of the singularity of the formally extrapolated running coupling 
to small momenta. However, the method of integration allows one to generalize the 
renormalon expressions to higher orders of the coupling evolution. In higher orders there 
is a freedom of choosing renormalization schemes for the coupling and one can always 
find a scheme in which the evolution trajectory is integrable in the IR region [^. With 
an integrable coupling one obtains a particular way of resummation of renormalon-type 
contributions. In QED the generalization to higher orders means that one inserts the 
full one-particle irreducible polarization function into the photon propagator. It is close 
in spirit to organization of the expansion in skeleton diagrams. In QCD, the renormalon 
chains beyond the leading order can be generated from the running coupling using Dyson- 
Schwinger equation p6| . 

Practical applications of the renormalon-type estimates are most frequent in the stan- 
dard QCD [^. An important consequence of the bubble-chain based analysis for heavy 
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quarks is an IR sensitivity of the pole mass |2^. Because the IR region is a nonpertur- 
bative one and the results based on the consideration of bubble chains are not rigorous 



e.g. [29, pO]) there is no unique quantitative evaluation of such a sensitivity. In high 



precision analysis of heavy quark properties it is generally assumed that the uncertainty 



of the heavy quark mass due to the IR sensitivity is numerically small I^I], |32 |. 

In the following we explore some possibilities to define an integral of the running 
coupling over the IR region by using a freedom of renormalization scheme (RS) choice 
in high-order PT. This provides a particular method of summation of asymptotic series 
related to the IR renormalons in QCD. For the method to work one should consider higher 
order PT corrections to the standard bubble-chain diagrams. 

First we consider an example of how higher order PT contributions can be used to 
give a definition of the IR renormalon contributions. Let the /3-function of a theory in 
some particular renormalization scheme be summed into the expression 

with the standard asymptotic /^-function of the form /?"*(q;) = — a^. The renormalization 
group equation for the running coupling constant a{z) 

z^a{z) = (3{a) (2) 

has a solution 



a{z) = — , = In — . (3) 

This solution has no singularities at positive z. However, the asymptotic solution 

has a pole at z = A^. Thus, the particular way of summing an infinite number of 
perturbative terms for the /^-function in eq. (|l|) has cured the Landau pole problem in 
a sense that the coupling becomes smooth in the IR region and the integration of the 
evolution trajectory can be performed explicitly. Let us stress that there are no nonPT 
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terms added but the freedom of choosing a renormahzation scheme for an infinite series 
defining the running couphng was used instead. This is, however, beyond the formal 
framework of perturbation theory where only finite order polynomials in the coupling 
constant are considered for /3-functions. The problem of RS ambiguity in high orders of 
PT was discussed in refs. |3^, ^ Q where several methods to reduce this ambiguity 
have been proposed. In contrast to previous approaches where the main criterion for a 
good scheme is a weak dependence on the change of the scheme parameterization, in the 
present paper we use the RS freedom to find a scheme that allows one to perform the 
integration over the infrared region without encountering a singularity. 

Eq. (|l]) can be considered either as a pure PT result in some particular renormahzation 
scheme after an infinite resummation or as an approximation of some exact /^-function (as 
obtained from a numerical study on the lattice, for instance) that might include nonpertur- 
bative terms as well. The essential point for us here is that the running coupling obeying 
the renormahzation group equation with such a /3-function can smoothly be extrapolated 
into the IR region. The running coupling given in eq. (§) has a correct asymptotic be- 
havior at 2; ^ oo. It is a regular expansion parameter for physical observables that obey 
the dispersion relation because it has no Landau pole on the physical cut. 

The generalization of eq. ([^) to QCD is now straightforward. Keeping in mind that 
first two coefficients /?o and f3i of the QCD /3-function are renormahzation scheme invariant 
one could suggest a possible "improved" /5-function as 

with the running coupling constant a{z) implicitly determined by the equation 

/^o^^T^ = ~ — R ( \ ~ kl3Qa[z) (6 
a{z) (3o \ (Jia{z) J 

The /3-function given in eq. (^) is bounded at large a and eq. (H) has a solution for a{z) 

that is defined on the positive semi-axis and has no singularities. 

The solution given in eq. (^) can be rewritten through an intermediate energy scale n 

without using the A parameter (the A parameter has no special meaning as a position of 



the pole of the running couphng). For the sake of simphcity we do that for the leading 
order expression given in eq. and find 

= , , (7) 

1 + aLq - ka^ + ^{1 + aLq - ka^ f + Aka"^ 

where Lq = ln((5^/ /i^) and a = a{^'^). This form of the evolution of the coupling constant 
is a particular regularization of the Landau singularity. This regularization introduces no 
imaginary part on the positive semi-axis (cf. ref. [0). Thus, in the pure PT framework 
one can eliminate the Landau pole of the coupling using the freedom of the RS choice 
in high orders. The smooth coupling requires no nonperturbative effects to regularize 
the integrals in which it appears. The integral over the small momenta determines a 
particular sum for a series related to an (extended) IR renormalon contribution. 

The expression given in eq. (|^) for the /5-function is not polynomial. Still even for a two- 
loop /5-function there can exist an IR fixed point that allows one to make an extrapolation 
to small energies. Indeed, for nj light flavors in QCD the first two coefficients of the (3- 
function read 

47r/3o = 11 - ^nj, {Ajifp^ = 102 - 



and for 



one has 



153 33 

<nf< 



19 ' 2 



/3o > 0, /3i < . 

Therefore, for n/ = 9, for instance, there exists an infrared fixed point and an extended 
IR-renormalon asymptotic series (generated with a two-loop /3-function) can be explicitly 
summed according to our prescription, i.e. in the second order approximation for the 
evolution of the coupling in QCD. 

The coupling constant defined in a particular renormalization scheme is not an immedi- 
ate physical quantity. However, one can write a formal definition of the effective coupling 
through some physical quantity. This definition can implicitly include nonperturbative 
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terms. For instance, taking the expression for Adler function 

2 f Re+e- {s)ds 



[s + Q 



2\2 



one can define 



and 38, 39 



1 



TT 



r^+e- is) = — arctan | - — — 



e(s) + e(s + A2)e(-s). 



(8) 



The unphysical singularity (a cut along a part of the negative semi-axis) corresponds to 
the unphysical pole in D{Q^). Subtracting this singularity in a minimal way we find the 
expression 



1 



(9) 



ln(Q7A2) Q2_A2 

that is regular for positive but has power corrections at large and is not polynomial 
in terms of the asymptotic charge 



rf'^^(Q2) = (ln(QVA2)) ' 



Indeed, for large one finds 



d^'iQ^ 



'(<5^) + . . . 



The corresponding /3-function written in terms of the asymptotic charge ci"** = ^ 

d 



(10) 



+ e-«(l-e~r 



also has an explicit nonperturbative term unlike eq. @. It is unclear whether the /3- 
function in eq. (|TTD can be rewritten in terms of the full charge (^(Q^) from eq. (P). 

Thus, there are various possibilities to eliminate the leading order pole and, thereby, to 
give an interpretation to the renormalon series as an integral of the running coupling over 
the IR region. A pure perturbative possibility means that higher order PT corrections 
make the evolution well-defined till the origin. This can be achieved with both finite-order 
and resummed expressions for the /3-function. Other possibilities can be nonperturbative 
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ones when the evolution is made smooth in the IR region by some nonperturbative pro- 
cedure. 

Within the IR-renormalon approach the power corrections emerge as a result of in- 
tegration of some sub-blocks of a whole diagram over the region of small momenta in 
the phase space of the diagram. The main generating sub-block in the IR-renormalon 
approximation is taken to be the running coupling. However, for defining the IR contri- 
butions there is no need to know the point-by-point evolution of the coupling in the IR 
domain. One has to define some integrals only. Finding only integrals is less demanding 
than determining the point-by-point evolution of the coupling in the IR domain. For 
phenomenological applications the information about a particular way of extrapolation of 
the coupling constant evolution into the IR domain can be encoded through just a nonPT 
parameter describing the leading order nonPT correction. In other words any particular 
prescription of summing the asymptotic series generated by the (extended) IR-renormalon 
chains corresponds to a numerical value of this parameter. For instance, the parameter 
k in eq. ([l|) describes a way of summing the IR-renormalon asymptotic series by using 
higher order corrections and changing the renormalization scheme in which the evolution 
is described. For phenomenological applications, however, it is convenient to retain a 
contact with the leading order asymptotic charge. By choosing the principal value (PV) 
prescription for the regularization of integrals of the asymptotic charge over the IR region 
one obtains a particular way of resumming the IR-renormalon contributions. This way 
is not unique. The residual freedom is parameterized by a nonPT quantity which is IR 
sensitive only The parameterization is written in the form 

Poa'^^is) = /3oa.(s) Ipv + AA^6is) + BA^6'is) + ... (12) 

The normalization is chosen such that for large s 

"-^^^ = /3oln(W) - ^''^ 
Eq. (|l^) is understood in a sense of distributions and as(s)|pv is the asymptotic charge 
containing the pole regularized by taking the principal value. Such a distribution is a 
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natural generalization of the classical solution to the renormalization group equation with 
the asymptotic /3-function P°-'^{x) = —PqX^. The generalized solution as{s)\py is not 
unique and allows some arbitrariness at the point s = that should be a series of 6- 
function and its derivatives. Quantities A,B, ... parameterize this arbitrariness. We have 
chosen the support of "nonperturbative terms" at s = instead of s = which suffices 
for the discussion in the present paper. However, if eq. ( |l2l) is integrated with a weight 
function that is singular in the origin one has to retain the exact location of the support 
s = A^ of the "nonperturbative terms". It is necessary, for instance, in case of finding 
renormalon corrections to the pole mass of a heavy quark. Note that multiplication is 
not defined for distributions therefore the non-leading PT corrections are supposed to be 
included in the leading order term in eq. (0) that results in the substitution A^ A^s 
(e.g. 1^1], ^ |4^)- The integration of expression in eq. (0) over the IR region with some 
numerical values for the parameters A and B should give the same result as an integration 
of the solution to eq. with some numerical value for the parameter k. For a while we 
retain two parameters in eq. (p!^). 

Now we consider several models for ^'^■^■^(s) and compute numerical values for param- 
eters A and B. In the spirit of the operator product expansion A = (there are no 
gauge invariant operators of dimensionality two) and B ~ (G^), where (G^) is the gluon 
condensate. With the constraint A = the parameters B is related to the parameter k 
from eq. (|^). Within such a model one expects that the pattern of extrapolation of the 
running coupling to the infrared region can be found by comparing the numerical value 
for the parameter B with that of the gluon condensate. The gluon condensate is con- 
nected with the violation of dilatation symmetry which is closely related to dimensional 
transmutation and appearance of the scale A. Note that quark condensates are related 
to chiral symmetry breaking and we do not consider them within our approach. 

Assuming a universality of extrapolation of the running coupling into the IR region we 
fix the parameter B from comparison with OPE. Taking the process of e+e~ annihilation 
where the operator product expansion is applicable we find the relation between the B 
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parameter and the gluon condensate. From the expression for Adler's function at large 
Euchdean 



one obtains the (naive) asymptotic spectral density in the leading order 



(14) 



Re+e-[s) = IH \ — — d [S] (15) 

TT iZ 



with 



°"'''-) = /?„ln(./A.V) - 

Note that 7r^-corrections are ignored (cf eq. (|])). Assuming the interpretation of the 
asymptotic charge a^"^ (s) in the above equations as a distribution on the semi-axis 
< s < oo with PV regularization we have 

<^^s) = af^-{sU + ^{9'sG')6'{s) . (17) 

This equation fixes a particular extrapolation of the coupling evolution into the IR domain. 
The chosen extrapolation is represented through the PV regularized asjTiiptotic charge 
and the gluon condensate contribution which has support only in the IR domain. This 
extrapolation is coordinated with OPE for the two-point correlator in e^e~ annihilation. 
Comparing eq. (|l^ with eq. (|12D one finds for three light quarks 

that translates into the following expression for the scale parameter Ag+e- through the 
gluon condensate 

For the standard numerical value {g^G'^) = 0.474 GeV^ we find 

546 

Ae+e- = — MeV. (18) 

154 



In e+e annihilation the parameter Ag+e- is an effective scale of the process [42, 41]. The 



use of the effective scale allows one to absorb the next-to-leading order PT correction into 
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al ^ (s). The numerical value of the effective scale is rather sensitive to the process and 
can characterize the scale where PT is violated and nonPT term are important numerically 



43|, W% - To find the relation between Ag+e- and the standard parameter h.^^ we consider 



the first perturbative correction to the coupling constant a^^'^ (s) in eq. (|I^) 

<rM.^SM+,7^!M)VoK») (19) 

where 

-iic(3) + „,(-H + !c(3) 

Here (^(-2) is Riemann (^-function. For the standard QCD with ri/ = 3 we find numerically 

W-"?'M+1.64&V + 0(aJ) (20) 



TT TT \ TT 



that leads to 

Ae+e- = 1.44Ams. (21) 
Using eqs. ( p!8| , |2T| ) one obtains a numerical prediction 

= 380 B~-^ MeV . (22) 

For the present experimental value of A^ extracted from the data on the r-lepton decay 
|5| 

A^ = 349 ± 61 MeV (23) 



B=[ y = 0.74 ^3.03. (24) 



one finds 

B={ 

.349 ±61. 

Below we consider several explicit models of extrapolation of the coupling evolution into 
the IR region that allow one to determine a numerical value for the parameter B. We find 
that the numerical value for the parameter B obtained from these models agrees with the 
result given in eq. (0). For the sake of technical simplicity we do not use sophisticated 
expressions for running coupling in the IR region stemming from RS redefinition. We 
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limit ourselves to models with rather a similar pattern of the IR behavior but expressible 
in simple analytical terms that makes the numerical analysis more transparent. 
Model 1. 

Switching the interaction off completely at small momenta: 
The system of equations for determining A and B is 



li{a) + A = 0, li{a^)-B = 0, 
where li{a) is a special function 



li{a) 



" dt 







ln(t) 

with the PV (principal value) prescription for the integration across the pole at t = 1 
for real positive a > 1 [0. The numerical solution to the system of equations with an 
additional constraint A = reads a = 1.45, B = /z(2.1) = 1.19. We retain only three 
significant figures.* 
Model 2. 

Freezing the running coupling at small momenta (e.g. p6| ) 
The system of equations for determining A and B is 

Qj (X^ 

li{a) + A, — |-— ^ = U{a^) - B. 



ln(a) ' ' 21n(a) 

The numerical solution with the constraint A = reads a = 3.85, B = 2.6.« 
Model 3. 

Regularization of the low energy behavior by the minimal subtraction of the pole singu- 
larity 



ln(s/A2) S-A2' 
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The system of equations to determine the parameters A and B has the form 

li{a) — ln(a — 1) = li{a) + A, li{a?) — a — ln(a — 1) = /i(a^) — B. 

The solution with A = is a = 2, 5 = 2. Note that this model has recently got some 
special attention |Q. • 
Model 4. 

Renormalization scheme modification of the /5-function and smooth extrapolation of the 
evolution to small momenta. The model is given by the pattern of running presented in 
eq. d^). The solution depends on the numerical value for the parameter k entering the 
expression for the /3-function in eq. (|I|). For k = 2 the numerical solution is a = 3.0, 
B = 3.0 which is close to that of Models 2,3. For larger k the solution is rather stable. 
For smaller values of k there may be no solution at all {k=l). This shows that the anzats 
chosen in eqs. ([T^JT7| ) restricts the possible extrapolations of the effective coupling into 
the IR region in Minkowskian domain and can be used to fix the numerical value for the 
B parameter through the gluon condensate.* 

Thus, the models of reasonably soft extrapolation of the running to small momenta give 



the value for the parameter B that agrees with the estimate given in eq. (@). Assuming 
that the numerical value for B has been found independently we can find a numerical value 
for the QCD scale parameter from eq. (0). Using the results of Models 1-3 B = 1.2 — 2.6 
or Bi = 1.05 - 1.27 one obtains = 300 ^ 362 MeV, or A^j^ = 331 ± 31 MeV that is 
in a reasonable agreement with the present data. In general, our result means that the 
numerical values of the gluon condensate and are compatible with each other for the 
smooth continuation of the running coupling into the IR region. The gluon condensate is 
interpreted as a quantity that corrects the PV-regularized leading-order running coupling 
at small momenta (in a sense of distributions) to provide a smooth continuation of the 
evolution into the IR domain. 

To conclude, we have exploited the RS freedom in high orders of PT to show that 
contributions of some infinite subsets of diagrams that are represented as integrals of 
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the running coupling over the IR region can exphcitly be summed in a generahzed way 
through the proper definition of the integraL The lack of any parameter or sensible strict 
criterion for choosing a particular set of diagrams that dominate high-order behavior 
beyond the leading order expression for the running coupling allows one to use the RS 
freedom for interpretation of the IR-renormalon singularities. If in a particular renor- 
malization scheme there are no singularities of the running coupling at small momenta, 
this scheme provides a particular recipe for the resummation of the asymptotic series 
related to the IR renormalons. As for phenomenological applications, the models with 
some extrapolation of the running coupling into the infrared region are used for practical 
calculation. Our results show that with an extrapolation chosen one cannot freely add 
the standard contribution of the gluon condensate to take into account nonperturbative 
effects and "improve" the computation - this contribution must be coordinated with the 
extrapolation of the coupling. 
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